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Abstract 

A major result about perturbations of integrable Hamiltonian sys- 
tems is the Nekhoroshev theorem, which gives exponential stability 
for all solutions provided the system is analytic and the integrable 
Hamiltonian is generic. In the particular but important case where 
the latter is quasi-convex, these exponential estimates have been gen- 
eralized by Marco and Sauzin if the Hamiltonian is Gevrey regular, 
using a method introduced by Lochak in the analytic case. In this pa- 
per, using the same approach, we investigate the situation where the 
Hamiltonian is assumed to be only finitely differentiable, for which it 
is known that exponential stability does not hold but nevertheless we 
prove estimates of polynomial stability. 

1 Introduction 

In this paper, we are concerned with the stability properties of near-integrable 
Hamiltonian systems of the form 

(H(e,i) = h(i) + f(6,i) 

\\f\<e«l 

where (0,1) E T n x M n are action- angle coordinates for the integrable part 
h and / is a small perturbation in some suitable topology defined by a norm 
| . | . More precisely, we are interested in the evolution of the action variables 
I(t), which are trivially constant in the absence of perturbation. 

1. The first main result in this direction is given by an application of the 
KAM theory on the persistence of quasi-periodic solutions (see |Pos01j for 
a recent exposition). Assuming h satisfies some non-degeneracy condition 
and the system is analytic, if e is sufficiently small there exists a constant c 
such that 

\I(t) - J 1 < cy/e, t G K, 
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for "most" initial actions Jo, more precisely for a set of large measure but 
with empty interior. When n = 2, this is even true for all solutions provided 
h is isoenergetically non degenerate, but for n > 3, the famous example of 
Arnold ( |Arn64j ) shows that there exist "unstable" solutions, along which 
the variation of the actions can be arbitrarily large no matter how small the 
perturbation is. From its very beginning, KAM theory was known to hold 
for non-analytic Hamiltonians (see |Mos62j in the context of twist maps). 
It is now well established in various regularity classes, including the C°° 
case (essentially by Herman, see |Bos86j and [Fej04j) and the Gevrey case 
( |Pop04| ). Following ideas of Moser, the theorem also holds if H is only of 
class C k , with k > 2n (see |Pos82j . [SalOi] . [5Z89] and also |Alb07j for a 
refinement), even though the minimal number of derivatives is still an open 
question, except in a special case for n = 2 ([He r86| ). 

2. Another fundamental result, which complements KAM theory, is given 
by Nekhoroshev's theorem ( |Nek77| . |Nek79| ) . If the integrable part h satis- 
fies some generic condition and the system is analytic, then for e sufficiently 
small there exist positive constants ci, 02,03, a and b such that 

\I{t) - J 1 < ci£ b , \t\ < c 2 exp(c 3 £- a ), 

for all initial actions Iq. Hence all solutions are stable, not for all time, 
but for an exponentially long time. In the special case where h is strictly 
quasi-convex, a completely new proof of these estimates was given by Lochak 
( |Loc92| ) using periodic averaging and simultaneous Diophantine approxi- 
mation. The method of Lochak has had many applications, in particular 
it was used by Marco and Sauzin to extend Nekhoroshev's theorem to the 
Gevrey regular case under the quasi-convexity assumption ([MS02 ). 

3. However, no such estimates have been studied when the Hamiltonian is 
merely finitely differentiable, and this is the content of the present paper. 
We will prove below (Theorem 12. ip that if H is of class C k , for k > 2, and 
h quasi-convex, then one has the stability estimates 

1 k — 2 

\I(t) — Jo I < C\£2^, \t\ < C2£~^"~ , 

for some positive constants c\ and C2, and provided that e is small enough. 
Of course, under our regularity assumption the exponential estimates have 
been replaced with polynomial estimates, and earlier examples show that 
exponential stability cannot possibly hold under such a weak regularity as- 
sumption (this is discussed in [MS04J. The proof will use once again the 
ideas of Lochak which, among other things, reduces the analytic part to its 
minimum and we will also follow the implementation of Marco and Sauzin 
in the Gevrey case. 

4. As we recalled above, KAM theory for finitely differentiable Hamiltonian 
systems has been widely studied, and so we believe that Nekhoroshev's esti- 
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mates under weaker regularity assumptions have their own interest. More- 
over, for obvious reasons, examples of unstable solutions (so-called Arnold 
diffusion) are more easily constructed in the non-analytic case, and it is 
a natural question to estimate the speed of instability (see [KL08aJ and 
|KL08b| for examples of class C k with a polynomial speed of diffusion). 
Finally, one of our motivations is to generalize these estimates using the 
method of |BN09j . where Lochak's ideas are extended to deal with analytic 
but C fc -generic unperturbed Hamiltonians, with k > 2n + 2. 



2 Main result 



1. Let T n = M n /Z n , and consider a Hamiltonian function H defined on the 
domain 

V R = T n X Br, 

where Br is the open ball of M n around the origin of radius R, with respect 
to the supremum norm | . | . As usual, we shall occasionally identify H with 
a function defined on M n x Br which is 1-periodic with respect to the first 
n variables. 

We assume that H is of class C k , for an integer k > 2, i.e. it is A;-times 
differentiable and all its derivatives up to order k extend continuously to the 
closure Vr. We denote by C k (T>R) the space of such functions, which is a 
Banach space with the norm 

\ H \c k (v E ) = sup SU P sup \d a H(x)\ 

0<l<k \a\=l \ xe V R J 

where x = (0, 1), a = (a>i, . . . , «2n) £ N 2n , \a\ = a\ + ■ ■ ■ + ct2 n and 

u — u l ... u 2n ■ 

In the case where the Hamiltonian H = h depends only on the action vari- 
ables, we will simply write \h\ C k^ BR y 

Our Hamiltonian H E C k (T>R) is assumed to be C fc -close to integrable, 
that is, of the form 

(H(0,I) = h(I) + f(8,I) 
\\f\c"{v H ) < e « 1 

where h is the integrable part and / a small perturbation of size e in the C k 
topology. 

We may assume that there exists a positive constant M such that 

\h\cHB R ) < M - (B) 
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Moreover, we will require that h satisfy the following quasi-convexity as- 
sumption: there exists a positive constant m such that 

V/GB fi ,V?)GR n , Vh(I).v = => V 2 h(I)v.v > m\v\ 2 . (C) 
2. Our main theorem is the following. 

Theorem 2.1. Let H be as in Jjj), with k > 3, h satisfying (B\) and < TCp . 

and sei 

k — 2 , 1 

a = , o = — . 

2n 2n 

T/ien i/iere exist eq, c\ and C2 such that if s < £o, a// solutions (6(t),I(t)) 
of H with 1(0) G -Br/2 satisfy 

\I(t)-I(0)\< Cl e b , \t\<c 2 e- a . 

First note that we have stated our theorem for of class C fc , > 3, but 
the statement also holds with no changes if H is of class C * -1 > 1 ) > 3, that 
is, is of class C k ~ 1 and its partial derivatives of order k — 1 are Lipschitz 
continuous. For C 2 systems, these estimates also hold true but it is not 
useful since trivially all solutions satisfy 

\m- i(o)i <s, \t\<i. 

In fact for C 2 Hamiltonians for which the derivatives up to order 2 are not 
more than continuous, we believe that one cannot obtain better estimates. 

Moreover, the time of stability obtained is "optimal" in the sense that 
one can construct examples of unstable orbits with a polynomial speed of 
diffusion, but we do not know what the optimal exponents should be. How- 
ever, using the geometric arguments of |BM10j one can easily improve the 
stability exponent a in order to obtain 

u _ o 

6, 



2(n-l) 



for 5 > but arbitrarily small. 

Let us finally point out that if H is C°°, then it is an immediate conse- 
quence of the above result that the action variables are stable for an interval 
of time which is longer than any prescribed power of e^ 1 , but even in this 
case exponential stability does not hold. 

3. As in the analytic or Gevrey case, we can also state a refined result near 
resonances. Suppose A is a sub-module of Z n of rank m, d = n — m and let 
Sa be the corresponding resonant manifold, that is 

S A = {I £ B R I k.Vh(I) = 0, k € A}. 

We can prove the following statement, which actually contains the previous 
one. 
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Theorem 2.2. Under the previous hypotheses, assume d(I(0), S&) < o~\fe 
for some constant a > 0, and set 

k-2 , 1 



Then there exist e' , c± and c' 2 such that if e < e' , one has 

\I(t)-I(0)\<c\e^, \t\ 

For A = {0}, d = n and S\ = Br/2, we recover Theorem 12,11 and 
therefore it will be enough to prove Theorem 12.21 

4. The constants eo, c\ and C2 depend only on h, more precisely they 
depend on k,n,R,M and m while the constants e' , c'i and c' 2 also depend 
on a and A. However we will not give explicit values for them in order 
to avoid complicated and rather meaningless expressions. Hence we shall 
replace them by the symbol • when it is convenient: for instance, we shall 
write u <• v when there exists a positive constant c depending only on the 
previous parameters, but not on /, such that u < cv. 

5. This paper is divided into two sections. The next section contains the 
analytical part of the proof, where we will construct a system of local coor- 
dinates for our Hamiltonian which is more convenient to study the evolution 
of the action variables. Then, in the last section we will conclude the proof 
using our convexity assumption and Dirichlet's theorem on simultaneous 
Diophantine approximation. 



3 Analytical part 

1. Given an action / 6 Br and denoting by oj = V/i(I) its frequency, we 
know from classical averaging theory that the relevant part of the perturba- 
tion 

f(9,I)= h(I)e l2wk£ 

is given by those harmonics associated with integers k G 7U 1 in resonance 
with oj, that is such that k.u = 0. Actually one can construct a symplectic, 
close-to-identity transformation $ defined around /, such that 

Ho$ = h + g + f 

where g contains only harmonics in resonance with u and / is a small remain- 
der. These are usually called resonant normal forms, and to obtain them 
one has to deal with small divisors k.oj which involve technical estimates. 
If the system is analytic, the above remainder / can be made exponentially 
small with respect to the inverse of the size of the perturbation, as was first 
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shown by Nekhoroshev. But for finitely difierentiable systems one might 
guess that the remainder can only be polynomially small, even though this 
should be difficult (or at least technical) to prove using the usual approach. 

2. It is a remarkable fact discovered by Lochak ( |Loc92| ) that to prove 
exponential estimates in the quasi-convex case with the analyticity assump- 
tion, it is enough to average along periodic frequencies, which are frequencies 
ui such that Tu G Z n \ {0} for some T > (see also [BN09] for an exten- 
sion of this method for generic integrable Hamiltonians) . These periodic 
frequencies correspond to periodic orbits of the unperturbed Hamiltonian, 
hence in this approach no small divisors arise. As a consequence this special 
resonant normal form is much easier to obtain. The aim of this section is 
to construct such a normal form, up to a polynomial remainder. This will 
be done in 13.31 But first we will recall some useful estimates concerning the 
C k norm in 13.11 and then prove an intermediate statement in 13.21 

3.1 Elementary estimates 

3. Let us begin by recalling some easy estimates. Given two functions 
f,g G C k (Vji), the product fg belongs to C k (Vn) and by the Leibniz rule 

\fd\c k (v R ) <• \f\c k (v R )\g\c k (v R )- 

The Poisson Bracket {f,g} belongs to C k ~ l (T>R), and by its definition and 
the Leibniz rule one gets 

\{f,9}\ C k ~ 1 (T> R ) <' \f\c k {V R )\9\c k (V R )- 

The above implicit constants depend only on n and k (in fact in the first 
estimate one can trivially modify the definition of the C k norm so as to have 
a constant equal to one, but this will not be important for us). These are 
very elementary facts, but we shall also need estimates concerning vector 
fields, canonical transformations and compositions. 

4. First, given a vector-valued function F G C k (T>R, W), F = (Fx, ...,F{) 
and I G N, we extend the norm component- wise, that is 

\ F \c k (V R ) = sup \Fi\ ck( p R y 

l<i<l 

Now for a function f G C k (T>R), we define its Hamiltonian vector field Xf 
by 

X f = (dif, -d e f) 

where 

d!f = (d h f, . . .,d In f), d e f = (d e J, . . .,d e J). 
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Obviously X f G C k ~ l (V R , R 2n ), and trivially 

\ x f\c k - 1 (v a ) < \f\c k {V R )- 

Moreover, by classical theorems on ordinary differential equations, if Xf is of 
class C k ~ 1 then so is the time-i map of the vector field Xf, when it exists. 
Assuming \def\c°{v R ) < r f° r some r < R (for example l-X^l^o^^) < r), 
then by the mean value theorem 

= X' i? _ r — >. V R 

is a well-defined C fc_1 -embedding. In the case where / is integrable, one can 
choose r = 0. 

In the sequel, we will need to estimate the C k norm of in terms of the 
C k norm of the vector field Xf. More precisely we need the rather natural 
fact that is C fc -close to the identity when Xf is C fc -close to zero. This 
is trivial for k = 0. In the general case, this follows by induction on k using 
on the one hand the relation 

= Id + f X f o 
Jo 

and on the other the formula of Faa di Bruno (see |AR67| for example), 
which gives bounds of the form 

\FoG\ ck <\F\ ck \G\ k ck 

and also 

iFoGb^-lF^lGl^ + lFl^lGl^ 

for C k vector-valued functions on appropriate domains (once again, the 
above implicit constants depend only on k). Let us state this as a lemma, 
for which we refer to [DH09J, Lemma 3.15 and appendix C, for a detailed 
proof. 

Lemma 3.1. Let X f G C k (V R , R 2n ), assume that \Xf\ c o^ R ^ < r and 

\ x f\c k (v R ) < !■ U) 

Then 

\<& f - Id\ C k(p R _ r ) <• \Xf\ C k(p R y 

The above implicit constant depends only on k and R. Now if H G 
C k (T>ji), under the above hypotheses we have H o G C k (T>R- T ) and the 
estimate 

\H o &\cnv R _ r) < \ h \chv r )\&\c* { -d r _ t) ( 2 ) 
follows trivially from the Faa di Bruno formula. 
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3.2 The linear case 



Following [MS02], we change for a moment our setting and we consider a 
perturbation of a linear Hamiltonian, more precisely the Hamiltonian 

H(e,i) = i(i) + f(ej) 

y** j 

\f\c"{v p ) < P « 1 

where p > is fixed and 1(1) = u.I is a linear Hamiltonian with a T-periodic 
frequency u. Recall that this means that 

T = inf{t >0\tuj £Z n \ {0}} 

is well-defined. In this context, our small parameter is \i. 

In the proposition below, we will construct a "global" normal form for 
the Hamiltonian which we will use in the next section to produce a 

"local" normal form around periodic orbits for our original Hamiltonian 

Proposition 3.2. Consider H as in it**]) with k > 2, and assume 

1><-1, |w|<l. (3) 
Then there exists a C 2 symplectic transformation 

$ : V p/2 -> V p 
with |$ — Id^m ) <-T/i suc/i i/iai 

^o(D = ; + 5 + /, 
with {g,l} = and i/ie estimates 

\g\cHv p/2 ) <• M> l/lc2(© p/2 ) < (Tp) k ~ 2 p 

hold true. 

First note that {g, 1} = means exactly that dgg.co = 0, and expanding g 
in a Fourier series, one easily sees that it contains only harmonics associated 
with integers k satisfying k.oo = 0. Therefore the above proposition gives 
indeed a resonant normal form, up to a polynomial remainder. 

Note also that we need our transformed Hamiltonian H o and hence 
our transformation <&, to be at least of class C 2 , simply because we need our 
transformed vector field to be of class C 1 to have existence and uniqueness 
of solutions (or course, a C 1 ' regularity would have been enough). This 
explains the factor (k — 2) in our stability exponent a. 
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Proof. First, as we are assuming \u>\ <• 1, we have 

< 1- ( 4 ) 

Our transformation $ will be obtained by a finite composition of averaging 
transformations. Let us define 

= P 
T 2(k - 2) 

and for j G {0, . . . ,k — 2}, let 

Pj = P~ jr > p/2. 

Then we claim that for any j G {0, . . . , k— 2}, there exists a C^ - - 7 symplectic 
transformation <3?j : — > D p with \$>j — ~[d\c k ~i(v p .) <:Tfi such that 

Ho$j =l + gj + fj, 

with gj and /j of class C fc_J , {<7j, /} = and the estimates 

\9j\c k -J(v Pj ) < M, \fj\c*-3(p p .) < {TpYp- 

The proposition follows easily by taking $ = <3?fc_2, 5 = gk-2 and / = /fc-2- 
We will prove the claim by induction on j G {0, . . . , k — 2}. For j = 0, 
there is nothing to prove since we can simply write H = I + go + /o with 
5o = 0, /o = / and therefore $o is the identity. Now assume the claim is 
true for some j G {0, . . . , k — 3}, and consider 

fT,- = Ho$ j = l + g j + f j . 

Let us define 

[/;] = -/ fjo^dt 

and 

*j = ^j[ '1.0 X>'"i>!</'- 

We have 

|[/j]lc*-j(x> Pj ) - l/j ° ^tlc*-^©^.) 
hence by ([2]) and @ we obtain 

I [/?] Ic fc -J {v Pj ) <• I fj I c* - i (D Pj ) 
and with our hypotheses of induction, this gives 

\[fj]\ck-J(v Pj ) < (Tp) 3 p. 
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Similarly 

\Xj\c^{V Pj )<T(T^ j H=(Tfi) j+1 

and since T[i <• 1, 

\Xj\c k -i{v p .) <Tfi. 

If we let $ Xj be the time-one map of the Hamiltonian vector field gen- 
erated by Xji then we will show that the map 

satisfies the assumptions. 

Indeed, thanks to the condition T/x <• 1, we can ensure that $* J , which 
is of class C is a well-defined embedding 

Moreover, as \X Xj \ C k-j-i(p \<-Tfj, and using once again T/i<-l, we can 
arrange condition (pQ) and apply Lemma 13.11 to obtain 

|*x* -Id| c . k - i - 1(J , />j+i) <r/z. 

Now 

" ld \c^-Hv P]+1 ) = l*j ° $ Xj - $ Xj + $ Xj - Id|c. k -i- 1(73pj+i) 

< l(*j-Id) * w lc*-i-i(23 PJ+1 ) 
+ |*^-Id| c *-j-i (17p . +i) 

<• l*j - Hlc'-i- 1 ^) + \® Xj - Id lc fc -^- 1 (© Pj+1 ) 

where we have used ([5]) in the last line. By our hypotheses of induction, this 
eventually gives 

- Id|o*-i-i(^ +1 ) < T ^ 
Now by Taylor's formula with integral remainder, we can expand 

H j+1 = H o = Z + + f j+1 

with 

ft+i = ft + ifil /i+i = / {ft + /]> Xj> ° 

J 

where /j = tfj + (1 — t)[fj], as one can check by a standard calculation. 
Since {gj, 1} = by our hypothesis of induction and obviously {[fj], 1} = 0, 
we have {gj+i, 1} = together with the estimate 

\9j+l\c k -J(V Pj ) - \9j\c*>-i(V p .) + \[fj]\c k -i(D Pj ) 

<• H 
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using Tfj, <• 1. For the remainder, using ([2]) and the trivial estimate for the 
Poisson bracket we compute 

\fj+l\ck-j-i(v Pj+1 ) < \{9j + fj, Xj} ° ®t } \c*-*- 1 (Pp j+1 ) 
<• \{9j + fpXj}\ck-i-i(v P .) 

< \{9j,Xj}\c k -i- 1 (V p .) + \{fj->Xj}\c k -i- 1 (V p .) 

<• \9j\c k -i{V P] )\Xj\c k -J(V P] ) + I /j I C fc -J (X> Pj . ) I Xj I C fe -J (X> Pj . ) 

using once again T\i <• 1. This concludes the proof. □ 
3.3 Normal form 

Now let us come back to our original setting which is the Hamiltonian 

(H(8,I) = h(I) + f(6,I) 

\\f\C(PB) < £ « L 

We will say that an action T 6 Br is T-periodic if its frequency vector 
V/i(I#) is T-periodic. In the proposition below, we will fix a T-periodic 
action T, I will be the linear integrable Hamiltonian 1(1) = oj.I associated 
with the periodic frequency uj = V/i(T) and we denote by 

11/ : TP xB R ^B R 

the projection onto the action space. 

Proposition 3.3 (Normal form). Suppose H is as in fl*]j, with h satisfy- 
ing IB\) . Under the previous hypotheses, let fx > be such that 

EK-/J 2 , /x<-l, T(i<l. (5) 

Then there exists a C 2 symplectic transformation 

$ : TP x B(h,Li) -> TP x B(h,2fi) 

with |ri/<I> — ldi\cofB(i t ,fj,)) <T/i 2 such that 

Ho^ = h + g + f, 

with {g,l} = and the estimates 

If + /Ic°(T«xB(I„m)) <• i^i l^fii/lc (T™xB(/»,At)) < V 2 
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Let us immediately explain how such coordinates will be used in the 
sequel. If we denote them by (9,1), then writing down the equations of 
motion for H = H o <£, using the fact that dgg.u = and the mean value 
theorem, one shows that I(t) remains close to the hyperplane orthogonal 
to Lj, passing through 1(0), for an interval of time governed by the size of 
dgf. Hence any potential drift has to occur along that hyperplane. Then 
this picture persists and gets only slightly distorted when we come back to 
the original coordinates (0,1), since the projection of <3? onto action space is 
close to identity. 

Proof. First note that since h satisfies ([B]), then 

\u\ = |V/t(I*)| < M 

and therefore the second part of condition ([3]) is satisfied. 

Without loss of generality, we may assume h(I*) = 0. To analyze our 
Hamiltonian H in a neighborhood of size fi around I*, we translate and 
rescale the action variables using the map 

tT li :(6,I»)^(6,I) = (e,I* + nI») 

which sends the domain T>2 = T n x B2 onto T n x B(I*, 2fi), and note that 
by the condition fi <• 1 , we can assume that the latter domain is included in 
V R . Let 

H„ = o<7„) 
be the rescaled Hamiltonian, so is defined on T>2 and reads 

H^O, J") = ^ 1 H(6, h + ^) = pT X h(h + fill*) + fT 1 f(e, h + 

for (9, II 1 ) £ T>2- Now using Taylor's formula we can expand h around /* to 
obtain 

h(h + fil 1 *) = fioj.Il* + fi 2 f (1 - t)V 2 h(h + tfil^I^.I^dt 

Jo 

= fiuj.l^* + /i 2 V( /M ) 

where we set 

h„(Ii*) = [ (1 - t)V 2 h(h + tfil^.I^dt. 
Jo 

Therefore we can write 
with 

U = ^ + ^(foa^). 
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Now we know that \f\c k (v R ) < £ <• /U 2 from the first part of condition © 
and \h\ C k(x, R ^ <■ 1 since we are assuming (fB]), so we obviously have 

\ffi\c k (V 2 ) < M- 

Therefore conditions ([3]) are satisfied and we can apply Proposition 13. 2| 
with p = 2, to the Hamiltonian = I + /„: there exists a C 2 symplectic 
transformation ^ : Di — > D2 with |<J>^ — Idj^^^ <• Tfj, such that 

# M o $ M = Z + + 
with {5^, /} = and the estimates 

Mc 2 (©i) <• M, l//Jc 2 (£>i) <• (7» fe ~V 
Moreover, if we set 

we have 1} = and |/ / i|c* 2 (X'i) <' A 4 ; an d so the transformed Hamiltonian 
can also be written as 

Hp o $ M = I + fihn + + / M . 

Now scaling back to our original coordinates, we define = o o cr^ 1 , 
therefore 

$ : TP x 5(7* — ► T n x B(J„2/z) 

and 

= n(l + + + / M ) o a" 1 

Observe that (/iZ + n 2 h^) o o" M 1 = /i, so we may set 

g = pff,oa-\ / = M / M oo-- 1 , 

and write 

iJo$ = Hj + /. 
It is obvious that {5, /} = with 

lfi , lc (T™x J B(/»,/i)) < M|/mIc°(I»i) 

and similarly 



13 



so 



Moreover, as dzf 



Iff + f\c°(T n xB{h,n)) < fJ 2 - 

lidgfn then 



dgf\c°(r n xB(h,fj,)) < C^VO 2 L j2 



and finally 



11/$ - Idj|cfO(fl(/, iA4 )) <-T/x 2 . 



is trivial. This ends the proof. 



□ 



4 Proof of Theorem [23] 



1. Now we can complete the proof of our Theorem 12.21 in the spirit of 
Lochak, following three elementary steps that use successively some arith- 
metic (simultaneous Diophantine approximation), some analysis (normal 
forms around periodic orbits) and some geometry (quasi-convexity). The 
analysis has been done in the previous section, and the arithmetic and ge- 
ometry are exactly the same as in the analytic case or the Gevrey case. 
Therefore instead of rewriting proofs which are well-known, we will merely 
explain the ideas and state the relevant results which can be found in |Loc92j 
and |MS02j . 

2. Let us begin with the arithmetic part, since in order to use our Proposi- 
tion [33] (the normal form), we will need to show that any action Iq G Br/2, 
which is close to some resonant surface 5a, can be approximated by a peri- 
odic action. Let uq = Vh(Io), then using the isoenergetic non-degeneracy of 
h, (which is easily implied by the quasi-convexity assumption), it is enough 
to approximate luq by a periodic vector oj. If A has rank n — 1, this is totally 
obvious, since necessarily we have A = oj 1 - n TP 1 for some periodic vector 
u! and so each action in S\ is periodic. Now in the case where A has rank 
m = n — d with d > 1, a good approximation is given by a theorem of 
Dirichlet, which moreover gives an explicit bound on the period T. 

Proposition 4.1. Let Iq £ Br^, A be a sub-module ofZ n ofrankm = n—d, 
with d > 1 and Q be a real number such that 



Q>1. 



(6) 



Then there exists a T-periodic action I* G Br such that 



Iq - I*| <• max d{I , S A ),T~ 1 Q ~ 



) 



and the period T satisfies 



1<T< Q. 
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The proof of the above proposition can be found in [MS02j, Corollary 

3.2. 

3. Now it remains to explain how the quasi-convexity hypothesis ([C]) on 
h, together with the normal form obtained in Proposition 13.31 wm enable 
us to control the variation \I(t) — 2o|> f° r an initial action Iq close to some 
periodic action I*. The idea goes as follows. 

Let u) = V/i(I*). We have already explained after Proposition 13.31 how in 
those new coordinates, the evolution of the actions in the direction generated 
by uj is controlled for some interval of time: if F is the hyperplane orthogonal 
to oj, then I(t) remains close to the affine subspace Iq + F for an interval of 
time |t| < t where t is essentially given by the inverse of the size of /. Now 
by preservation of energy, for all time I(t) remains close to the unperturbed 
energy hypersurface Eq = {I £ Br \ h(I) = h(Io)}, and as the latter is 
strictly convex by quasi-convexity of h, the connected component of Iq in 
Eq n (Jo + F) is in fact bounded and so is the variation \I(t) — Iq\ for \t\ < r. 

This idea is formalized in the proposition below. Once again, is a 
T-periodic action and I is the linear Hamiltonian with periodic frequency 
u = Vh(F). 

Proposition 4.2. Under the previous hypotheses, let r > 0, r > and 

H = h + g + f£ C 2 (T n x B(F,r)) 
with h satisfying IU\) , {g, 1} = and the estimates 

\g + f\c°{T n xB(h,r)) < r 2 , \d§f\c°(T*xB(I*,r)) < r 2 T~ X . 

If 

r <■ 1 (7) 
then for any initial condition (9o,Io) S T n x B(I*,r), the solution satisfies 

\I(t) - I | <-r, \t\ < r. 

Once again, we refer to [MS02) . Corollary 3.1, for a complete proof. 

4. Let us now conclude the proof of Theorem 12.21 In a first step we will use 
Proposition 14.11 to find a periodic action close to our initial action, then in 
a second step we will apply Proposition 13.31 to find adapted coordinates and 
the third step will consist in applying Proposition 14.21 to control the evolu- 
tion of the action variables in those coordinates, and hence in the original 
coordinates. 

Proof of Theorem \2.2[ Let Iq £ Br^, let A be a sub-module of 7L n of rank 
m = n — d, and assume that d(Io, S\) <■ Ve (note that this last assumption 
is void if A is trivial). 
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First step. In the case d = 1, any action J € £a has a frequency V/i(J) 
which is a multiple of some non zero vector l; t £ Z", therefore we can choose 
a periodic action J* £ Sa so 

1 7 - J* I = d{I ,S A )< ye 

and the period T trivially satisfies T<1. In the case d > 1, we apply 
Proposition 14.11 with 

d-l 

Q = -e 2d . 

and the condition ([6]) gives a first smallness condition on e. Observe that 
L_ j_ 

Q d-i = . £ 2d } hence the periodic action J* given by the proposition satisfies 

| J -/*|<- max (d(I ,S A ),T- 1 Q~^ <T~ 1 e^ 
and the period T 

KT<e~Tr. (8) 

Second step. Having found a periodic action, we will now apply Propo- 
sition 13.31 with 

With this choice, for d > 1 the first part of condition ([5]) is satisfied thanks 
to the upper bound (JSj) on the period T (for d = 1, this is trivial). The other 
conditions in ([5]) give only further smallness conditions on e. Applying the 
proposition, we have a C 2 symplectic transformation 

$ : T n x B{h,n) -»■ T n x B(J», 2/i) 

with |II/<3? — Id|co(B(/» ft)) <" J"^ 2 such that 

H o§ = h + g + f 

with {g, /} = and the estimates 

Is + /lc°(T»xB(J,,Ai)) < /A l^/lco^x-BC/*,^)) <• {Tv) k ~ 2 H 2 - 

Let us write H = H o & and (0, J) are the new coordinates in T n x J? (J*, //). 

Third step. Now we set 

r = - M , r = -(T^) fc - 2 , 

and we apply Proposition 14.21 to the Hamiltonian H. To do so, we need 
to impose condition ([7]) and this gives our last smallness condition on e. 
Therefore we find 

|J(i)-/o|<M, \t\<(T^) k - 2 , 
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and recalling that 



[1 = -T l £Zd <• £2d 

this gives 

— — 1 k — 2 

\I(t) - I \ <£2d, |t|<£"23-. 

Now since 

|H/$ - Id| c o (B(/ , iAt)) <• Tii 2 <• , 
standard arguments give the conclusion 

1 k — 2 

- io| <£2d, |i| <- e~2d-. 

This ends the proof. □ 

Acknowledgements. The author thanks Jean-Pierre Marco, Laurent Nie- 
derman and also H. Scott Dumas for a careful reading of a first version of 
the manuscript. 
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